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Updated values of running quark and lepton masses at GUT (Grand unified theories) scales are 
important for fermion mass model building, and to calculate neutrino masses, in GUTs . We present 
their values at GUT scales, in SM, MSSM and 2HDM theories, using the latest values of running 
quark and lepton masses. 
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■ I. INTRODUCTION 

o : 

Gauge theories are very attractive theories to explain the origin of all interactions among the fundamental particles. 
Standard model (SM) is a gauge theory based on group SU(2)l x U(l)y x SU(3)c (G213). In SM, all fundamental 
particles get their masses via the celebrated Higgs mechanism. One of the major goals of current research in exper- 
imental and theoretical high energy physics is to understand the origin of all fermion masses and mixings, including 
those of neutrinos. Although SM has been very successful in explaining many of the observed experimenatl results, 
some questions remain unanswered in it. Gauge hierarchy problem, unification of gauge couplings, neutrino masses, 
origin of baryon asymmetry of the Universe (BAU), being the most importane ones. Some of these problems can be 
circumvented if we consider two higgs doublet model (2HDM), minimal supersymmetric standard model (MSSM) [1], 
Q_i and GUTs. Although problem of gauge hierarchy is not solved by 2HDM, unification of gauge couplings is possible 
at GUT scales, in MSSM, and also after embedding them in non-SUSY GUTs like SO(10) 0. Very recently 0, we 
, have shown unification of the three gauge couplings a\y (for U(1)y), c*2L (for SU(2)l), and u^c (f° r SU(3i)c) in 
non-SUSY SM with additional flavor symmetries, and also estimated limits on proton life time. 

It is now a welll established fact that neutrinos have mass, and mix with each other and oscillate to other flavors. We 
know that in SM, neutrinos masses can not be explained, and hence we need to go to theories beyond standard model 
(BSM). One of the most promising theories, to explain small neutrino masses, is the grand unified theory (GUT), like 

^\ • SO(10), in which all the fermions, inlcuding the right handed (RH) neutrino, are present in a single I6-dimensional 
representation. These theories require running masses and mixings of quarks and charged leptons at GUT scales, for 

0^ , calculating neutrino masses. In theories based upon quark-lepton unification, like L-R symmetric SU(2)l x SU(2)r x 
SU(4)c group, these values are also required at intermdeiate scales. Unification of fundamental forces is based upon 

\Q , gauge symmetries which contain the standard model with fermions in the fundamental representations. Thus, the 

■ explaination of fermion masses and mixings must emerge from a successful unified gauge theory. And hence, the 
running fermion masses are required to build underlying textures and models for existence of appropriate unified 

I ■ theory. 

J> ' Values of running masses of quarks and charged leptons at higher scales in SM, 2HDM and MSSM are available in 
, literature 0. They have been used quite extensively, by many researchers, e.g. in 

• 0) for constructing neutrino masses 

• 0j f° r studying structures of unified thoeries 

• 0to study type II seesaw dominance in Non-SUSY and split SUSY SO (10) theory 

• 0, for study of SO(10) models, to explain fermion masses and mixing angles, including neutrino masses. 

• Q3 for study of inverse seesaw in NonSUSY SO(10) theories 

But, in all these works, older values from (J|have been used, and new data for fermion masses are available, for using 
as input at lower scales. The aim of present work is, to update these values, and fill the gap. We have used latest 
data for masses and couplings from PDG [ll|. Conversion of MS to DR scheme is done using formulas given in 
(l2^ . and top quark mass is taken from [r|- Following the anlysis of 0, we use RGEs for Yukawa couplings, gauge 
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couplings and VEVs separately, and calculate running values of fermion masses at GUT scale. These values at other 
intermediate scales, calculation of neutrino masses using them, will be presented elsewhere [14j . 

The paper has been organized as follows. In Section 2, we give a pedagogical discussion on fermion masses. Section 
3 contains methodology of, how to run fermion masses from one energy scale to another. In Section 3, our new results, 
on updated values of running fermion masses, at GUT scale 2 x 10 16 GeV, have been presented. Discussions and 
conclusions have been given in Section 5. 

II. A PEDAGOGICAL DISCUSSION ON FERMION MASSES 

Now, we will have a pedagogical discussion on fermion masses. According to quantum field theory (QFT), the 
"bare" masses in the Lagrangian are infinite for all particles, but divergent loop contributions to the propagator 
cancel them out to give finite "dressed" masses. This is called renormalization. These dressed particle masses are 
actually measured in experiments. So in the case of an electron, for example, the experimentally measured electron 
mass is an input parameter to the theory, and according to QFT, the bare electron mass must be infinite, but the 
mass "runs" from infinity at very small length scales, to a constant at very large length scales ("IR fixed point"). So 
this IR-limit value is the same as the experimentally measured value. 

We know that quarks are confined, and free quarks cannot be observed experimentally. This short distance con- 
finement is believed to be because of nonperturbative effects, and is associated with the scale Aqcd ~ 2 GeV 2 . At 
energies greater than Aqcd, the QCD is perturbative. Since free quarks do not exist at energy scales less than Aqcd 
(also called infrared (IR) limit) , mass for them is not well defined. Hence quark masses are scale dependent, and they 
are aften defined at a energy scale. The scale dependent quark masses are called 'current' or 'running' quark mass, 
and they are renormalization scheme dependent. But equivalence of these renormalization scheme-dependent quark 
masses can be established with renormalization group equations (RGEs). The 'constituent' quark mass is believed 
to be roughly the mass that contributes to observed mass of hadron, for example. Nonrelativistic quark models use 
constituent quark masses, the constituent mass of up and down quarks are ~ 350 MeV. 

For quarks masses also 'running' takes place, but instead of converging to a constant, they diverge at the energy scale 
Aqcd- They become infinite at a much smaller length scale. This makes perfect sense because quarks are confined 
into hadrons and can't be observed macroscopically.The masses given in PDG [ll|are the values of the 'running' 
masses at some energy scale greater than (length scale smaller than) Aqcd, defined in some specific renormalization 
scheme. 



III. RUNNING OF MASSES AND COUPLINGS USING RGES 

In the renormalization theories, where the Yukawa couplings and the VEVs run separately [l5[-[23j , the Dirac mass 
of a fermion can be defined as 

M i (jj) = Y i (j i )v i (p). (1) 

Here, Mi(fj,) is the Dirac mass of the i-type fermion, Yi(/i) is corresponding Yukawa coupling, and Vi(fi) is the running 
VEV (Vaccum expectation value), at the scale [i. In these scenarios, the Yukawa couplings and VEVs run separately, 
independent of each other. Many authors have used these [15]-[23], see [4] (Das, Parida) for a complete discussion. 
The relevant terms of the Lagrangian, for masses of fermions, in SM, can be written as: 

L = qiXucpUR + qiY D (pd R + l L Y E <t>e R + h.c. (2) 

Here, (f> is the higgs particle, v(fi) its running VEV at scale fi, is the left handed quark doublet, ur is the right 
handed quark, d R is the right handed quark, is the left handed lepton doublet, and cr is the right handed electron. 
Since in SM , no right handed neutrinos are present, there is no term in the Lagrangian for the neutrino mass. Similarly, 
for 2HDM and MSSM, this can be written as: 



Here, 



L = qhYu4>uUB, + qLY D (f> D dR + IlYe^dcr + h.c. 
< (j) v >= vu(fj.) = v((j,)sin/3, < 4* D >= Vd(/*) = v(fi)cos/3 



(3) 



(4) 
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and 

tan/3((j,) = vu(p)/vd(p) (5) 

Now, we write the RGEs for running of Yukawa and gauge couplings, for the SM, 2HDM, and MSSM, along with 
their RG coefficients. They have been given in [J], but we present them here for the sake of c omp leteness only. The 
one-loop RGEs for Yukawa couplings, for SM, MSSM and 2HDM, can be written as [H|-(I1], HiMH 

16tt 2 ^ = [Tr(3YuY^ + 3aY D Y^ + aY E Y^) 



dt 



hbYuY^ + cY D Yl) £ '<#]Y V (6) 



dY 

16tt 2 — £ = [TriiaYuYl + WnYl + YEYl) 



dt 



\(bY D Yl + cYjjY^) -Y,C\ d) i\Y D (7) 



16^ 2 ^ = [Tr{^aY v Y^ + 2Y D Yl+Y E Yl) 



dt 



\bY E Yl-Y^C\ e) 9 ^]Y E (8) 



The RGEs for the VEV in SM, upto 2-loop have been derived using wave-function renormalisation of the scalar field 
[15-16, 18-19, 21-22 ], and the 1-loop equation is 

16^* = c (v) g 2 _ Tr{? Y vY \ + WdY 1 + Ye yI)]v (9) 

i 

Here, t — ln/i. The RGEs for v a (a = u, d) in the 2HDM up to 1-loop and MSSM up to 2-loops are available in [15-18, 
20]. The 1-loop equations in both theories are 

16tt 2 ^ = [£» C^gf - Tr(3YuY^v u (1Q) 



16^2^ = [E . Cfg 2 _ Tr{3YD Yl + Y E Yl)v u (n) 

The RGE for the gauge couplings for the three models are 

W^=b i9 f (12) 

2-loop contributions are available in literature [ 15-18, 21-26], and we use them from Das, ParidaQ. 

Using above RGEs, we run the values of fermion masses, from low scale Mz to higher scale 2 x 10 16 GeV. The input 
values of running fermion masses at Mz have been taken from PDG (ill ], and (l2l ]. Our results have been presented 
in next section. 

IV. RESULTS 

The new results of our computations have been presented in Tables (I- VI ). We have presented comparisons of all 
our results with older values (Das, Parida, EPCJ 2001). We have used mass of the Higgs to be 125 GeV. It can be 
noted that from a recent global analysis [23], mass of the Higgs boson has been expected to be around this value. 
The scale of supersymmetry breaking, Ms = 1 Tev has been used. It is worth mentioning here that some signatures 
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Table I: Comparison of input values at Mz 



Fermion 


Mass (This work) 


Nlass (Das, Parida) 


m u 


1 /I £+0.56 T\/T^Af 

l-45_ 45 MeV 


A.66_ a 45 MeV 


m c 


cot; _|_ oa AJfcAT 


0/ i el iviev 


mt 


1/z.l ± U.D ± U.y OeV 


iqi _i_iQ 
ioi ± io Oev 


TTld 


z.yu_Q ^ iviev 


4 - b9 -0.66 MeV 


m 3 


MeV 


93.4+^;® MeV 


m b 


2.82lg;S GeV 


3.0 ± 0.11 GeV 


m e 


0.486847 MeV 


0.486847 MeV 


«V 


102.75138 ± 0.00033 MeV 


102.75138 ± 0.00033 MeV 


m T 


1.74669t°;°°£? GeV 


1.74669±g:°8S? GeV 



of SUSY have been detected at LHC in third family of fermions (2^|. The pole mass of top quark is used from PDG 
[U, to be m t = 172.9 ± 0.6 ± 0.9 GeV. This is first converted to running mass m t (M z ) = 172.1 ± 0.6 ± 0.9, as 
described in Xing et al |l2j. This value is used for SM and 2HDM. Then, for MSSM only, we convert this running 
value m t (Mz) to DR(dimensional regularization) scheme value, by using Eq. (22) of Xing et. al. [12j, and find this 
to be m t {Mz)DR = 169.9 ± 0.6 ± 0.9. The latest PDG value l/a(M z ) = 128.91 and a s (M z ) = 0.1189 ± 0.0020 are 
used in our analysis. 



A. Running fermion masses in SM at GUT scale = 2 x 10 lf> GeV 
TABLE-II COMPARISON OF FERMION MASSES IN SM, 2-LOOP 



Fermion 


Update Mass (This work) 


Mass (Das, Parida) 


m u 


0.4565±8;133| MeV 


0.8351±g;i^g MeV 


m c 


0.2225l[™ GeV 


0.2426±8;g|^ GeV 


m t 


70.5188 + g-|2?| GeV 


75.43481^^ GeV 


m d 


1.0773±g;ti^ MeV 


1.7372lg;«« MeV 


m s 


20.4323i|:2i2 MeV 


34.5971±f;f|fl MeV 


m b 


0.9321±g;{j$Jj GeV 


0.9574±°;°°jg GeV 


m e 


0.4413 ± 0.0003 MeV 


0.4413 ±0.0001 MeV 


m M 


93.116 T 0.0117 MeV 


93.143lt{J;{£25 MeV 


m T 


1.6109 T 0.00003 GeV 


1.5834t°;°°° 5 GeV 



B. Running fermion masses in MSSM at GUT scale = 2 x 10 16 GeV 
TABLE-III COMPRISON OF MASSES IN MSSM, 2-LOOP, TAN^= 10 



fermion 


mass (this work) 


mass (Das, Parida) 


m u 


0.3961±8;i|g| MeV 


0.7238+°^ MeV 


m c 


0-19301™^ GeV 


0.21031°:™° GeV 


m t 


T1.0883±i;|H§ GeV 


82.4333l^;^6 GeV 


m d 


0.9316l°;^s MeV 


1.5036l°;g35 MeV 


m s 


17.6702tte95o MeV 


29.94541^5444 MeV 


nib 


0.9898±8;gl§^ GeV 


1.06361°;^ GeV 


m e 


0.3585l°;°™2 MeV 


0.3585 ± 0.0003 MeV 


rrifj, 


75.639 T 0.0003 MeV 


75.6715l8:o50i McV 


m T 


1.314fi±g:88S GeV 


1.2922±°;°™ GeV 
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TABLE-IV COMPARISON OF MASSES IN MSSM, 2-LOOP, TAN£=55 



fermion 


mass (this work) 


mass (Das, Parida) 


m u 


0.3963+°, ^ MeV 


0.7244+°'^^ MeV 


m c 


0.1932+° °240 GeV 


0.2105+° °^} GeV 

— U.Uzll 


mt 


80.4472+^H GeV 


95.1486+^-^° GeV 


m-d 


0.9284 + °'SS MeV 


1.4967lg;Sre MeV 


m s 


17.6097+4;^ MeV 


29.8135+J™ MeV 


nib 


1.24241°;°^° GeV 


1.4167t°+*° 3 4 GeV 


m e 


0.3569;°;°°°i MeV 


0.3565;°;°°°^ MeV 




75.35707.°;°™ MeV 


75.2938^^515 MeV 




1.6459t°;°206 GeV 


1.6292t°;°^ GeV 



C. Running fermion masses in 2HDM at GUT scale = 2 x 10 b GeV 

TABLE-V COMPARISON OF MASSES IN 2HDM, 1-LOOP, TAN^=10 



fermion 


mass (this work) 


mass (Das, Parida) 


m u 


0.47761 MeV 


0.874918 im MeV 


m c 


0.2328±g;g|gg GeV 


0.2542l8;8 2 255 GeV 


m t 


74.1053l^°^ 3 GeV 


79.6373l;! 1 1 ^ 74 GeV 


m d 


1.12741°;^ MeV 


1.820418:2743 MeV 


m s 


21.3821+° °°°° MeV 


36.2544l^;°™8 MeV 


mb 


1.1615™ GeV 


1.2309l8;8?38 GeV 


m e 


0.4407t°;°°° 3 3 MeV 


0.4407;°;°°°1 MeV 


rrifj, 


92.9898;o;oiig MeV 


93.0197;8.oi22 2 MeV 


m T 


1.61277^:88^ GeV 


1.585ll8;° 88 5 5 GeV 


COMPARISON OF MASSES IN 2HDM, 1-LOOP 


fermion 


mass (this work) 


mass (Das, Parida) 


m u 


0.4776l8;^9 6 MeV 


0.874918: \™l MeV 


m c 


0.2328+°;°^ GeV 


0.254218:8^ GeV 


m t 


77.3752±J;£g GeV 


83.9317±£j™ GeV 


ma 


1.1274l° j°^ MeV 


1.820418:^3 MeV 


m s 


21.3836t^8s3 McV 


36.25841*;°™ MeV 


m b 


1.3053t°;°f 5 ° GeV 


1. 412818: U62 GeV 


m e 


0.4407t° °°° 3 3 MeV 


0.4407;8:8ooi MeV 


rrifj. 


93.0222^8'8m MeV 


93.0536^ oi36 MeV 


m T 


1-8138+«;S GeV 


1.785ll8:8i8 6 7 GeV 



', TAN/3=55 



V. DISCUSSIONS AND CONCLUSIONS 

We have presented updated values of running fermion masses in SM, 2HDM and MSSM at GUT scale, at tan/3 = 10 
and 55, using 2-loop RGEs for SM and MSSM, and 1-loop RGEs for the 2HDM. It can be seen from our results (Tables 
I- VI ) that these new values of fermion masses are quiet different from their older counterparts (Das, Parida Q). 
They can be used for calculation of neutrino masses in GUTs at higher scales, as well as for buliding of theories for 
fermion mass models. Here, we would like to mention that we have verified our calculations, by reproducing the values 
reported in Das, Parida [4]. Also, our values are different from values reported in Xing. et. al. [l^. This is beacuse 
we have used a different scheme for running fermion masses from low scale to GUT scale. We have used RGEs for 
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Yukawa couplings and VEVs separately. As discussed in text, fermion masses using this scheme Qhave been used 
extensively in literature. Hence the results presented in this paper are very important. The values of running fermion 
masses at other intermediate scales, and calculation of neutrino masses using them will be presented elsewhere 
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